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Devices for nano- and molecular size electronics are currently a focus of re- 
search aimed at an efficient current rectification and switchingQ A few generic 
molecular scale devices are reviewed here on the basis of first-principles and 
model approaches. Current rectification by (ballistic) molecular quantum dots 
can produce the rectification ratio < 100. Current switching due to confor- 
mational changes in the molecules is slow, on the order of a few kHz. Fast 
switching (~1 THz) may be achieved, at least in principle, in a degenerate 
molecular quantum dot with strong coupling of electrons with vibrational ex- 
citations. We show that the mean-field approach fails to properly describe 
intrinsic molecular switching and present an exact solution to the problem. 
Defects in molecular films result in spurious peaks in conductance, apparent 
negative differential resistance, and may also lead to unusual temperature and 
bias dependence of current. The observed switching in many cases is extrinsic, 
caused by changes in molecule-electrode geometry, molecule reconfiguration, 
metallic filament formation through, and/or changing amount of disorder in a 
molecular film. We give experimental examples of telegraph "switching" and 
"hot spot" formation in the molecular films. 



1 Introduction 

Current interest in molecular electronics is largely driven by expectations 
that molecules can be used as nanoelectronics components able to comple- 
ment/replace standard silicon CMOS technologyfTJ [2] on the way down to 
<10nm circuit components. The first speculations about molecular electronic 
devices (diodes, rectifiers) were apparently made in mid-1970s 3_. That origi- 
nal suggestion of a molecular rectifier has generated a large interest in the field 
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and a flurry of suggestions of various molecular electronics components, es- 
pecially coupled with premature estimates that silicon-based technology can- 
not scale to below l/xm feature size. The Aviram-Ratner's Donor-insulator- 
Acceptor construct TTF— a— TCNQ (D+ — a — A~ ^ see details below), where 
carriers were supposed to tunnel asymmetrically in two directions through in- 
sulating saturated molecular cr— 'bridge', has never materialized, in spite of ex- 
tensive experimental effort over a few decades ^4) . End result in some cases ap- 
pears to be a slightly electrically anisotropic insulator, rather than a diode, un- 
suitable as a replacement for silicon devices. This comes about because in order 
to assemble a reasonable quality monolayer of these molecules in Langmuir- 
Blodgett trough (avoiding defects that will short the device after electrode 
deposition) one needs to attach a long 'tail' molecule C18 [= (CH2)i8] that 
can produce enough of a Van-der-Waals force to keep molecules together, but 
C18 is a wide-band insulator with a bandgap Eg ^ 9 — lOeV. The outcome 
of these studies may have been anticipated, but if one were able to assemble 
the Aviram-Ratner molecules without the tail, they could not rectify anyway. 
Indeed, a recent ab-initio study[5] of D^(tA~ prospective molecule showed 
no appreciable asymmetry of its I-V curve. The molecule was envisaged by 
EUenbogen and Love [6] as a 4-phenyl ring Tour wire with dimethylene insulat- 
ing bridge in the middle directly connected to Au electrodes via thiol groups. 
Donor-acceptor asymmetry was produced by side NH^ and NO^ moieties, 
which is a frequent motif in molecular devices using the Tour wires. The rea- 
son of poor rectification is simple: the bridge is too short, it is a transparent 
piece of one-dimensional insulator, whereas the applied field is three dimen- 
sional and it cannot be screened efficiently with an appreciable voltage drop 
on the insulating group in this geometry. Although there is only 0.7eV en- 
ergy separation between levels on the D and A groups, one needs about 4eV 
bias to align them and get a relatively small current because total resonant 
transparency is practically impossible to achieve. Remember, that the model 
calculation implied an ideal coupling to electrodes, which is impossible in 
reality and which is known to dramatically change the current through the 
molecule (see below) . We shall discuss below some possible alternatives to this 
approach. 

It is worth noting that studies of energy and electron transport in molecu- 
lar crystals started already in early 1960s. It was established in mid-1960s in 
what circumstances charge transport in biological molecules involves electron 
tunneling [8]. It was realized in mid-1970s that since the organic molecules are 
'soft', energy transport along linear biological molecules, proteins, etc. may 
proceed by low energy nonlinear collective excitations, like Davydov solitons 
[9] (see review [lO]). 

To take over from current silicon CMOS technology, the molecular elec- 
tronics should provide smaller, more reliable, functional components that can 
be produced and assembled concurrently and are compatible with CMOS 
for integration. The small size of units that molecules may hopefully pro- 
vide is quite obvious. However, meeting other requirements seems to be a 
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very long shot. To beat alternative technologies for e.g. dense (and cheap) 
memories, one should aim at a few Tb/in^ (> 10^^ — 10^^ bit/cm^), which 
corresponds to linear bit (footprint) sizes of 3 — 10 nanometers, and an oper- 
ation lifetime of ~ 10 years. The latter requirement is very difficult to meet 
with organic molecules that tend to oxidize and decompose, especially under 
conditions of very high applied electric field (given the operational voltage 
bias of ~ IV for molecules integrated with CMOS and their small sizes on 
the order of a few nanometers). In terms of areal density, one should com- 
pare this with rapidly developing technologies like ferroelectric random access 
(FERAM)[TT] or phase-change memories (PCM) [12]. The current smallest 
commercial nano-ferroelectrics are about 400x400 nm^ and 20-150 nm thick 
[13], and the 128 x 128 arrays of switching ferroelectric pixels bits have been 
already demonstrated with a bit size < 50nm (with density ~Tb/in^)[14]. The 
phase-change memories based on chalcolgenides GeSbTe (GST) seem to scale 
even better than the ferroelectrics. As we see, the mainstream technology for 
random-access memory approaches molecular size very rapidly. At the same 
time, the so-called "nanopore" molecular devices [35] have comparable sizes 
and yet to demonstrate a repeatable behavior (for reasons explained below). 

In terms of parallel fabrication of molecular devices, one is looking at self- 
assembly techniques (see, e.g. [151 [H], and references therein). Frequently, the 
Langmuir-Blodgett technique is used for self-assembly of molecules on water, 
where molecules are prepared to have hydrophilic "head" and hydrophobic 
"tail" to make the assembly possible, see e.g. Refs. [T7l[T8]. The allowances 
for a corresponding assembly, especially of hybrid structures (molecules inte- 
grated on silicon CMOS), are on the order of a fraction of an Angstrom, so 
actually a picotechnology is required [2] . Since it is problematic to reach such 
a precision any time soon, the all-in-one molecule approach was advocated, 
meaning that a fully functional computing unit should be synthesized as a 
single supermolecular unit [2]. The hope is that perhaps directed self-assembly 
will help to accomplish building such a unit, but self-assembly on a large scale 
is impossible without defects [121 dS], since the entropic factors work against 
it. Above some small defect concentration ("percolation") threshold the map- 
ping of even a simple algorithm on such a self-assembled network becomes 
impossible[l9]. 

There is also a big question about electron transport in such a device con- 
sisting of large organic molecules. Even in high-quality pentacene (P5) crys- 
tals, perhaps the best materials for thin film transistors, the mobility is a mere 
1-2 cm^/V-s (see e.g. [20]), as a result of carrier trapping by interaction with a 
lattice and necessity to hop between P5 molecules. The situation with carrier 
transport through long molecules (> 2 — 3nm) is, of course, substantially dif- 
ferent from the transport through short rigid molecules that have been envis- 
aged as possible electronics components. Indeed, in short molecules the domi- 
nant mode of electron transport would be resonant tunneling through electri- 
cally active molecular orbital(s) [21j . which, depending on the workfunction 
of the electrode, affinity of the molecule, and symmetry of coupling between 
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molecule and electrode may be one of the lowest unoccupied molecular orbitals 
(LUMO) or highest occupied molecular orbitals (HOMO) [HI US] . Indeed, it 
is well known that in longer wires containing more than about 30-40 atomic 
sites, the tunneling time is comparable to or larger than the characteristic 
phonon times, so that the polaron (and/or bipolaron) can be formed inside 
the molecular wire |24j . There is a wide range of molecular bulk conductors 
with (bi)polaronic carriers. The formation of polarons (and charged solitons) 
in polyacetylene (PA) was discussed a long time ago theoretically in Refs. [25] 
and formation of bipolarons (bound states of two polarons) in Ref. [55] . Po- 
larons in PA were detected optically in Ref. [27] and since then studied in great 
detail. There is an exceeding amount of evidence of the polaron and bipolaron 
formation in conjugated polymers such as polyphenylene, polypyrrole, poly- 
thiophenc, polyphenylene sulfide [28], Cs-doped biphenyl n-doped bithio- 
phene polyphenylenevinylene(PPV)-based light emitting diodes [3T], and 
other molecular systems. Both intrinsic and extrinsic parameters play a role 
in determining the electrical and optical properties of polymer films: spatial 
range of 7r-electron delocalization, interchain interaction, morphology, amount 
of defects and disorder, carrier density, etc (for a brief review of carrier trans- 
port mechanisms and materials see Ref. [32]). 




The latest wave of interest in molecular electronics is mostly related to 
recent studies of carrier transport in synthesized linear conjugated molecu- 
lar wires (Tour wires [l]) with apparent non-linear I-V characteristics [neg- 
ative differential resistance (NDR)] and "memory" effects[33l jSU |35], vari- 
ous molecules with a mobile microcycle that is able to move back and forth 
between metastable conformations in solution (molecular shuttles) |36| and 
demonstrate some sort of "switching" between relatively stable resistive states 
when sandwiched between electrodes in a solid state device[37] (see also [38]). 
There are also various photochromic molecules that may change conforma- 
tion ("switch") upon absorption of light [39], which may be of interest to 
some photonics applications but not for the general purpose electronics. One 
of the most serious problems with using this kind of molecules is power dissi- 
pation. Indeed, the studied organic molecules rule, very resistive (in 
the range of ~ 1M]7— lGi7, or more). Since usually the switching bias exceeds 
0.5V the dissipated power density would be in excess of 10 kW/cm^, which is 
orders of magnitude higher than the presently manageable level. One can drop 
the density of switching devices, but this would undermine a main advantage 
of using molecular size elements. This is a common problem that CMOS faces 
too, but organic molecules do not seem to offer a tangible advantage yet. There 
are other outstanding problems, like understanding an actual switching mech- 
anism, which seems to be rather molecule-independent [38] , stability, scaling, 
etc. It is not likely, therefore, that molecules will displace silicon technology, 
or become a large part of a hybrid technology in a foreseeable future. 

First major moletronic applications would most likely come in the area of 
chemical and biological sensors. One of the current solutions in this area is 
to use the functionalized nanowires. When a target analyte molecule attaches 
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from the environment to such a nanowire, it changes the electrostatic potential 
"seen" by the carriers in the nanowire. Since the conductance of the nanowire 
device is small, even one chemisorbed molecule could make a detectable change 
of a conductance |40] . Semiconducting nanowires can be grown from seed metal 
nanoparticlesp^S] . or it can be carbon nanotubes (CNT), which are studied 
extensively due to their relatively simple structure and some unique properties 
like very high conductance j4lj. 

In this paper we shall address various generic problems related to electron 
transport through molecular devices, and describe some specific molecular 
systems that may be interesting for applications as rectifiers and switches, 
and some pertaining physical problems. We shall first consider systems where 
an elastic tunneling is dominant, and interaction with vibrational excitations 
on the molecules only renormalizes some parameters describing tunneling. 
We shall also describe a situation where the coupling of carriers to molecu- 
lar vibrons is strong. In this case the tunneling is substantially inelastic and, 
moreover, it may result in current hysteresis when the electron- vibron interac- 
tion is so strong that it overcomes Coulomb repulsion of carriers on a central 
narrow-band/conjugated unit of the molecule separated from electrodes by 
wide band gap saturated molecular groups like (CIl2)„, which we shall call 
a molecular quantum dot (molQD). Another very important problem is to 
understand the nature and the role of imperfections in organic thin films. It 
is addressed in the last section of the paper. 

2 Role of molecule-electrode contact: extrinsic molecular 
switching due to molecule tilting 

We have predicted some while ago that there should be a strong dependence 
of the current through conjugated molecules (like the Tour wires P]) on the 
geometry of molecule-electrode contact [22l [23]. The apparent "telegraph" 
switching observed in STM single-molecule probes of the three-ring Tour 
molecules, inserted into a SAM of non-conducting shorter alkanes, has been 
attributed to this effect [34]. The theory predicts very strong dependence of 
the current through the molecule on the tilting angle between a backbone of a 
molecule and a normal to the electrode surface. Other explanations, like rota- 
tion of the middle ring, charging of the molecule, or effects of the moieties on 
the middle ring, do not hold. In particular, switching of the molecules without 
any NO2 or NII2 moieties have been practically the same as with them. 

The simple argument in favor of the "tilting" mechanism of the conduc- 
tance lies in a large anisotropy of the molecule-electrode coupling through 
TT-conjugated molecular orbitals (MOs). In general, we expect the overlap and 
the full conductance to be maximal when the lobes of the p-orbital of the end 
atom at the molecule are oriented perpendicular to the surface, and smaller 
otherwise, as dictated by the symmetry. The overlap integrals of a p-orbital 
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Fig. 1. Schematic representation of the benzene-dithiolate molecule on top and 
hollow sites. End sulfur atoms are bonded to one and three surface gold atoms, 
respectively, 6 is the tilting angle. 

with orbitals of other types differ by a factor about 3 to 4 for the two ori- 
entations. Since the conductance is proportional to the square of the matrix 
element, which contains a product of two metal-molecule hopping integrals, 
the total conductance variation with overall geometry may therefore reach 
two orders of magnitude, and in special cases be even larger. 

In order to illustrate the geometric effect on current we have considered 
a simple two-site model with p— orbitals on both sites, coupled to electrodes 
with s— orbitals [23j . For non-zero bias the transmission probability has the 
resonant form ([5|) with line widths for hopping to the left (right) lead 
{Ffi) . The current has the approximate form (with F = Fl + Fn) 



where is the tilting angle O [22], Fig. [T] 

The tilting angle has a large effect on the I-V curves of benzene-dithiolate 
(BDT) molecules, especially when the molecule is anchored to the Au elec- 
trode in the top position, Fig.[2l By changing 9 from 5° to just 15°, one drives 
the I-V characteristic from the one with a gap of about 2V to the ohmic one 
with a large relative change of conductance. Even changing 9 from 10° to 
15° changes the conductance by about an order of magnitude. The I-V curve 
for the hollow site remains ohmic for tilting angles up to 75° with moderate 
changes of conductance. Therefore, if the molecule in measurements snaps 
from the top to the hollow position and back, it will lead to an apparent 




; sin 9, gV^ ^ -Elumo ^ -E-homOi 
cx sin^ 9, qV > -Elumo - £^homo, 



(1) 
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switching [53]. It has recently been reaUzed that the geometry of a contact 
strongly affects coherent spin transfer between molecularly bridged quantum 
dots [42j . It is worth noting that another frequently observed extrinsic mech- 
anism of "switching" in organic layers is due to electrode material diffusing 
into the layer and forming metallic "filaments" (see below). 




Fig. 2. Effect of tilting on I-V curve of the BDT molecule. Fig. [T] Current is in 
units of Iq = 77.5/xA, 6 is the tilting angle. 



3 Molecular quantum dot rectifiers 

Aviram and Ratner speculated about a rectifying molecule containing donor 
(D) and acceptor (A) groups separated by a saturated a-bridge (insulator) 
group, where the (inelastic) electron transfer will be more favorable from A to 
D [5]. The molecular rectifiers actually synthesized, C16II33 — 7Q3CNQ, were 
of somewhat different D — n — A type, i.e. the "bridge" group was conjugated 
[4]. Although the molecule did show rectification (with considerable hystere- 
sis), it performed rather like an anisotropic insulator with tiny currents on 
the order of f O^^^A/molecule, because of the large alkane "tail" needed for 
LB assembly. It was recently realized that in this molecule the resonance does 
not come from the alignment of the HOMO and LUMO, since they cannot be 
decoupled through the conjugated tt— bridge, but rather due to an asymmetric 
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voltage drop across the molecule [33]. Rectifying behavior in other classes of 
molecules is likely due to asymmetric contact with the electrodes [331 HS] , or 
an asymmetry of the molecule itself [3^. To make rectifiers, one should avoid 
using molecules with long insulating groups, and we have suggested using rel- 
atively short molecules with " anchor" end groups for their self-assembly on a 
metallic electrodes, with a phenyl ring as a central conjugated part [47]. This 
idea has been tested in Ref. [48] with a phenyl and thiophene rings attached 
to a (CH2)i5 tail by a CO group. The observed rectification ratio was <10, 
with some samples showing the ratio of about 37. 

We have recently studied a more promising rectifier like --S-(CH2)2-Naph— 
(CH2)io— S— with a theoretical rectification < 100[49], Fig.[3| This system has 
been synthesized and studied experimentally [50j . To obtain an accurate de- 
scription of transport in this case, we employ an ab-initio non-equilibrium 
Green's function method [5T] . The present calculation takes into account only 
elastic tunneling processes. Inelastic processes may substantially modify the 
results in the case of strong interaction of the electrons with molecular vibra- 
tions, see Ref. [52] and below. There are indications in the literature that the 
carrier might be trapped in a polaron state in saturated molecules somewhat 
longer than those we consider in the present paper [53j . One of the barriers 
in the present rectifiers is short and relatively transparent, so there will be no 
appreciable Coulomb blockade effects. The structure of the present molecular 
rectifier is shown in Fig. [Sj The molecule consists of a central conjugated part 
(naphthalene) isolated from the electrodes by two insulating aliphatic chains 
(CH2)„ with lengths Li(L2) for the left (right) chain. 




Fig. 3. Stick figure representing the naphthalene conjugated central unit separated 
from the left(right) electrode by saturated (wide band gap) alkane groups with length 

il(2)- 

The principle of molecular rectification by a molecular quantum dot is 
illustrated in Fig.|31 where the electrically "active" molecular orbital, localized 
on the middle conjugated part, is the LUMO, which lies at an energy A above 
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the electrode Fermi level at zero bias. The position of the LUMO is determined 
by the work function of the metal qcj) and the affinity of the molecule qx, 
A = Z\lumo — Qi4'~ x) ■ The position of the HOMO is given by zinoMO = 
Ai,\jMO—Eg, where Eg is the HOMO-LUMO gap. If this orbital is considerably 
closer to the electrode Fermi level Ep, then it will be brought into resonance 
with Ep prior to other orbitals. It is easy to estimate the forward and reverse 
bias voltages, assuming that the voltage mainly drops on the insulating parts 
of the molecule, 

VF = iil + 0, = + (2) 

Vf/Vr = ^ = L,/L2, (3) 

where q is the elementary charge. A significant difference between forward 
and reverse currents should be observed in the voltage range Vp < \V\ < Vr. 
The current is obtained from the Landauer formula 

I^^JdE ifiE) - J{E + qV)] giE, V). (4) 

We can make qualitative estimates in the resonant tunneling model, with the 
conductance g{E,V) = T{E,V)/q, where T{E,V) is the transmission given 
by the Breit-Wigncr formula 

T{E, V) = (5) 

Emo is the energy of the molecular orbital. The width /L(fl) ^ / D = 
/oe~^'^''^i(^' , where t is the overlap integral between the MO and the electrode, 
D is the electron band width in the electrodes, k the inverse decay length of 
the resonant MO into the barrier. The current above the resonant threshold 
is 

Ti 

We see that increasing the spatial asymmetry of the molecule {L2/L1) changes 
the operating voltage range linearly, but it also brings about an exponential 
decrease in current [17]. This severely limits the ability to optimize the rec- 
tification ratio while simultaneously keeping the resistance at a reasonable 
value. To calculate the I-V curves, we use an ab-initio approach that combines 
the Keldysh non-equilibrium Green's function (NEGF) with pseudopotential- 
based real space density functional theory (DFT) [FT]. The main advantages 
of our approach are (i) a proper treatment of the open boundary condition; 
(ii) a fully atomistic treatment of the electrodes and (iii) a self-consistent cal- 
culation of the non-equilibrium charge density using NEGF. The transport 
Green's function is found from the Dyson equation 

{G^y' = {GSy'~v, (7) 
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Fig. 4. Transmission coefficient versus energy E for rectifiers — S-(CH2)2- 
CioHe — (CH2)n— S— , n = 2,4,6,10. A indicates the distance of the closest MO 
to the electrode Fermi energy (Ef = 0). 



where the unperturbed retarded Green's function is defined in operator form 
as (G^) = {E + iO)S — H, H is the Hamiltonian matrix for the scatterer 
(molecule plus screening part of the electrodes). S is the overlap matrix, Si_j = 
iXilXj) for non-orthogonal basis set orbitals Xij and the coupling of the scat- 
terer to the leads is given by the Hamiltonian matrix V = dia,g[Si_i, 0, Sr.r], 
where I (r) stands for left (right) electrode. The self-energy part which 
is used to construct the non-equilibrium electron density in the scattering re- 
gion, is found from = — 2ilm [f{E)Si^i + f{E + qV)Er^r] , where £'/,i(r,r) 
is the self-energy of the left (right) electrode, calculated for the semi-infinite 
leads using an iterative technique [HT]. S'^ accounts for the steady charge 
"flowing in" from the electrodes. The transmission probability is given by 

TiE, V) = 4Tr [(ImEij) Gf;, (ImE,,.) G^j] , (8) 

where G^(^) are the retarded (advanced) Green's function, and E the self- 
energy part connecting left (/) and right (r) electrodes IFT, and the cur- 
rent is obtained from Eq. The calculated transmission coefficient T{E) 
is shown for a series of rectifiers — S-(CH2)m-CioH6— (CH2)„— S— for m — 2 
and n = 2, 4, 6, 10 at zero bias voltage in Fig. [4l We see that the LUMO is 
the molecular orbital transparent to electron transport, lies above Ep by an 
amount A = 1.2- 1.5eV. The transmission through the HOMO and HOMO-1 
states, localized on the terminating sulfur atoms, is negligible, but the HOMO- 
2 state conducts very well. The HOMO-2 defines the threshold reverse voltage 
Vr, thus limiting the operating voltage range. Our assumption, that the volt- 
age drop is proportional to the lengths of the alkane groups on both sides, is 
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quantified by the calculated potential ramp. It is close to a linear slope along 
the (CH2)„ chains [IHj- The forward voltage corresponds to the crossing of the 
LUMO(F) and hr{V), which happens at about 2V. Although the LUMO de- 
fines the forward threshold voltages in all molecules studied here, the reverse 
voltage is defined by the HOMO-2 for "right" barriers (CH2)n with n — 6, 10. 
The TV curves are plotted in Fig. [5l We see that the rectification ratio for 




Fig. 5. I-V curves for naphthalene rectifiers — S-(CH2)2-CioH6 — (CH2)n — S — , n — 
2, 4, 6, 10. The short-dash-dot curve corresponds to a cyano-doped (added group 
-C=N) n = 10 rectifier. 



current in the operation window /+ //_ reaches a maximum value of 35 for 
the "2-10" molecule (m = 2, n = 10). Series of molecules with a central single 
phenyl ring [47j do not show any significant rectification. One can manipu- 
late the system in order to increase the energy asymmetry of the conducting 
orbitals (reduce A). To shift the LUMO towards Ep, one can attach an elec- 
tron withdrawing group, like — C=N [49J. The molecular rectification ratio is 
not great by any means, but one should bear in mind that this is a device 
necessarily operating in a ballistic quantum-mechanical regime because of the 
small size. This is very different from present Si devices with carriers diffusing 
through the system. As silicon devices become smaller, however, the same 
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effects will eventually take over, and tend to diminish the rectification ratio, 
in addition to effects of finite temperature and disorder in the system. 

4 Molecular switches 

There are various molecular systems that exhibit some kind of current "switch- 
ing" behavior [331 [33 [3H]) "negative differential resistance" [55] . and 
"memory" [35] . The switching systems are basically driven between two states 
with considerably different resistances. This behavior is not really sensitive to 
a particular molecular structure, since this type of bistability is observed in 
complex rotaxane-like molecules as well as in very simple alkane chains (CH2),i 
assembled into LB films .38) . and is not even exclusive to the organic films. 
The data strongly indicates that the switching has an extrinsic origin, and is 
related either to bistability of molecule-electrode orientation [22l O [34] , or 
transport assisted by defects in the film [Ml [55] . 

4.1 Extrinsic switching in organic molecular films: role of defects 
and molecular reconfigurations 

Evidently, large defects can be formed in organic thin films as a result of 
electromigration in very strong field, as was observed long ago[5B]. It was 
concluded some decades ago that the conduction through absorbed [53] and 
Langmuir-Blodgett [55] monolayers of fatty acids (CH2)„, which we denote as 
Cn, is associated with defects. In particular, Polymeropoulos and Sagiv stud- 
ied a variety of absorbed monolayers from C7 to C23 on AI/AI2O3 substrates 
and found that the exponential dependence on the length of the molecular 
chains is only observed below the liquid nitrogen temperature of 77K, and 
no discernible length dependence was observed at higher temperatures [53]. 
The temperature dependence of current was strong, and was attributed to 
transport assisted by some defects. The current also varied strongly with the 
temperature in Ref. [55] for LB films on AI/AI2O3 substrates in He atmo- 
sphere, which is not compatible with elastic tunneling. Since the He atmo- 
sphere was believed to hinder the AI2O3 growth, and yet the resistance of the 
films increased about 100-fold over 45 days, the conclusion was made that the 
"defects" somehow anneal out with time. Two types of switching have been 
observed in 3-30 /im thick films of polydimethylsiloxane (PDMS), one as a 
standard dielectric breakdown with electrode material "jet evaporation" into 
the film with subsequent Joule melting of metallic filament under bias of about 
IQOV, and a low-voltage (<1 V) "ultraswitching" that has a clear "telegraph" 
character and resulted in intermittent switching into a much more conductive 
state [57]. The exact nature of this switching also remains unclear, but there is 
a strong expectations that the formation of metallic filaments that may even 
be in a ballistic regime of transport, may be relevant to the phenomenon. 
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Fig. 6. Experimental setup for mapping local conductance. AFM produces local 
deformation of top electrode and underlying organic film. The total conductance of 
the device is measured and mapped. (Courtesy C.N. Lau). 



Recently, a direct evidence was obtained of the formation of "hot spots" 
in the LB films that may be related to the filament growth through the film 
imaged with the use of AFM current mapping [58]. The system investigated in 
this work has been Pt/stearic acid (C18)/Ti (Pt/C18/Ti) crossbar molecular 
structure, consisting of planar Ft and Ti electrodes sandwiching a monolayer 
of 2.6-nm-long stearic acid (C18H360H) molecules with typical zero-bias re- 
sistance in excess of 10^ J7. The devices has been switched reversibly and re- 
peatedly to higher ("on") or lower ("off") conductance states by applying 
sufficiently large bias voltage Vt to top Ti electrode with regards to Ft coun- 
terelectrode, Fig. [SI 

Interestingly, reversible switching was not observed in symmetric Ft/C18 /Ft 
devices. The local conductance maps of the Ft/C18/Ti structure have been 
constructed by using an AFM tip and simultaneously measuring the current 
through the molecular junction biased to VJ, = O.IV (AFM tip was not used 
as an electrode, only to apply local pressure at the surface). The study re- 
vealed that the film showed pronounced switching between electrically very 
distinct states, with zero-bias conductances 0.17/iS ("off" state) and 1.45/iS 
("on" state). Fig. [Tj 

At every switching "on" there appeared a local conductance peak on the 
map with a typical diameter ~ 40nm, which then disappeared upon switching 
"off" , Fig. [7] (top inset) . The switching has been attributed to local conducting 
filament formation due to electromigration processes. It remains unclear how 
exactly the filaments dissolve under opposite bias voltage, why they tend to 
appear in new places after each switching, and why conductance in some 
cases strongly depends on temperature. It is clear, however, that switching in 
such a simple molecule without any redox centers, mobile groups, or charge 
reception centers should be extrinsic. Interestingly, very similar "switching" 
between two resistive states has also been observed for tunneling through thin 
inorganic perovskite oxide films [59 . 
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Voltage (V) 

Fig. 7. I-V characteristic showing the reversible switching cycle of the device 
(bottom inset) with organic film. The arrows indicate sweep direction. A negative 
bias switches the device to a high conductance state, while the positive one switches 
it to low conductance state. The mapping according the schematic in Fig. |6] shows 
the appearance of "hot spots" after switching (top inset). (Courtesy C.N. Lau). 

There have been plenty of reports on non-linear I-V characteristics hke 
negative differential resistance (NDR) and random switching recently for 
molecules assembled on metal electrodes (gold) and silicon. Reports on NDR 
for molecules with metal contacts ( Au, Hg) have been made in [331 EH HO] ■ 
It became very clear though that most of these observations are related to 
molecular reconfigurations and bond breaking and making, rather to any in- 
trinsic mechanism, like redox states, speculated about in the original Ref. 33j. 
Thus, the NDR in Tour wires was related to molecular reconfiguration with 
respect to metallic electrode[551 [33], NDR in ferrocene-tethered alkyl mono- 
layers [60] was found to be related to oxygen damage at high voltage[6T]. 
Structural changes and bond breaking have been found to result in NDR in 
experiments with STM [BH[S31[M] and mercury droplet contacts |65). 

Several molecules, like styrene, have been studied on degenerate Si surface 
and showed an NDR behavior [5S]. However, those results have been carefully 
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checked later and it was found that the styrene molecules do not exhibit NDR, 
but rather sporadically switch between states with different current while held 
at the same bias voltage (the blinking effect) [57]. 

The STM map of the styrene molecules (indicated by arrows) on the Silicon 
(100) surface shows that the molecules are blinking, see Fig. [51 The blinking 
is absent at clean Si areas, dark (D) and bright (C) defects. This may indicate 
a dynamic process occurring during the imaging. Comparing the panel (a) 
and (b) one may see that some molecules are actually decomposing. The 
height versus voltage spectra over particular points are shown in Fig. [51:. The 
featureless curve 1 was taken over a clean silicon dimer. The other spectra were 
recorded over individual molecules. Each of these spectra have many sudden 
decreases and increases in current as if the molecules are changing between 
different states during the measurement causing a change in current and a 
response of the feedback control, resulting in a change in height so there exists 
one or more configurations that lead to measurement of a different height. 
Evidently, these changes have the same origin as the blinking of molecules 
in STM images. Figure [Sb reveals clear structural changes associated with 
those particular spectroscopic changes. In each case where a dramatic change 
in spectroscopy occurred, the molecules in the image have changed from a 
bright feature to a dark spot. This is interpreted as a decomposition of the 
molecules. A detailed look at each decomposed styrene molecule, at locations 
3, 4, 5, and 6, shows that the dark spot is not in precise registry with the 
original bright feature, indicating that the decomposition product involves 
reaction with an adjacent dimer [67j. 

The fact that the structural changes and related NDR behavior are not 
associated with any resonant tunneling through the molecular levels or redox 
processes, but are perhaps related to inelastic electron scattering or other 
extrinsic processes, becomes evident from current versus time records shown 
in Fig. [9l Ref . [67] . The records show either no change of the current with time 
(1), or one or a few random jumps between certain current states (telegraph 
noise). The observed changes in current at a fixed voltage obviously cannot 
be explained by shifting and aligning of molecular levels, as was suggested in 
Ref. [68], they must be related to an adsorbate molecule structural changes 
with time. Therefore, the explanation by Datta et al. that the resonant level 
alignment is responsible for NDR does not apply [68]. As mentioned above, 
similar telegraph switching and NDR has been observed in Tour wires |34| 
and other molecules. Therefore, the observed negative differential resistance 
apparently has similar origin in disparate molecules adsorbed on different 
substrates, and has to do with molecular reconformation/reconfiguration on 
the surface. 



16 A.M. Bratkovsky 




VcltBg«{V) 

Fig. 8. STM images or styrene molecules on clean Si(lOO) before and after spec- 
troscopy over the area 75 x240 A. (Top left) Bias -2 V, current 0.7 nA. Only the 
styrene molecules (indicated by arrows) are blinking during imaging. The clean Si 
surface, bright defects (marked C), and dark sports (marked D) do not experience 
blinking. (Bottom left) bias -2 V, current 0.3 nA. STM image of individual styrene 
molecules indicated with numbers 2-6. Styrene molecules 3-6 have decomposed. De- 
composition involves the changing of the styrene molecule from a bright feature to 
a dark depression and also involves the reaction with an adjacent dimer. Styrene 
molecule 2 does not decompose and images as usual with no change of position. 
(Right) Height-voltage spectra taken over clean silicon (1) and styrene molecules 
(2-6). The spectra taken over molecules show several spikes in height related to 
blinking in the images. In spectra 3-6, an abrupt and permanent change in height 
is recorded and is correlated with decomposition, as seen in the bottom left image. 
Spectrum 2 has no permanent height change, and the molecule does not decompose. 
(Courtesy J. Fitters and R. Wolkow). 

4.2 Intrinsic polarization and extrinsic conductance switching in 
molecular ferroelectric PVDF 

The only well established, to the best of our knowledge, intrinsic molecu- 
lar switching (of polarization, not current) under bias voltage was observed 
in molecular ferroelectric block co-polymcrs polyvinylidene[69| . Ferroelectric 
polymer films have been prepared with the 70% vinylidene fluoride copoly- 
mer, P(VDF-TrFE 70:30), formed by horizontal LB deposition on aluminum- 
coated glass substrates with evaporated aluminum top electrodes. The poly- 
mer chains contain random sequence of (CH2)n(CF2)m blocks, fluorine site 
carries a strong negative charge, and in the ferroelectric phase most of carbon- 
fluorine bonds point in one direction. The fluorine groups can be rotated and 
aligned in very strong electric field, ~ 5MV/cm. As a result, the whole molecu- 
lar chain orders, and in this way the macroscopic polarization can be switched 
between the opposite states. The switching process is extremely slow, how- 
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Fig. 9. Variation of current through styrene molecules on Si(lOO) with time. Tun- 
neling conditions were set at -3 V and 0.05 nA. Abrupt increases and decreases 
in current relate to changes of the molecule during the spectroscopy. Some experi- 
ments show no changes in current (curve 1), others show various kinds of telegraph 
switching. (Courtesy J. Fitters and R. Wolkow). 

ever, and takes 1-10 seconds (!) '7U^,^7T. This is not surprising, given strong 
Coulomb interaction between charged groups and the metal electrodes, pin- 
ning by surface roughness, and steric hindrance to rotation. This behavior 
should be suggestive of other switching systems based on one of few monolay- 
ers of molecules, and other nontrivial behavior involved [71). 

The switching of current was also observed in films of PVDF 30 monolay- 
ers thick. The conductance of the film was following the observed hysteresis 
loop for the polarization, ranging from ~ 1 x 10^^ — 2 x 10~^J7~^[72]. The 
phenomenon of conductance switching has these important features: (i) It is 
connected with the bulk polarization switching; (ii) there is a large ~ 1000 : 1 
contrast between the ON and OFF states; (iii) the ON state is obtained only 
when the bulk polarization is switched in the positive direction; (iv) the con- 
ductance switching is much faster than the bulk polarization switching. The 
conductance switches ON only after the 6s delay, after the bulk polarization 
switching is nearly complete, presumably when the last layer switches into 
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alignment with the others, while the conductance switches OFF without a 
noticeable delay after the application of reverse bias as even one layer re- 
verses (this may create a barrier to charge transfer). The slow ~ 2s time 
constant for polarization switching is probably nucleation limited as has been 
observed in high-quality bulk films with low nucleation site densities[73]. The 
duration of the conductance switching transition ~ 2ms may be limited only 
by the much faster switching time of individual layers. 

The origin of conductance switching by 3 orders of magnitude is not 
clear. It may indeed be related to a changing amount of disorder for tunnel- 
ing/hopping electrons. It is conceivable that the carriers are strongly trapped 
in polaron states inside PVDF and find optimal path for hopping in the mate- 
rial, which is incompletely switched. This is an interesting topic that certainly 
in need of further experimental and theoretical study. 

5 Molecular quantum dot switching 
5.1 Electrically addressable molecules 

For many applications one needs an intrinsic molecular "switch" , i.e. a 
bistable voltage-addressable molecular system with very different resistances 
in the two states that can be accessed very quickly. There is a trade-off be- 
tween the stability of a molecular state and the ability to switch the molecule 
between two states with an external perturbation (we discuss an electric 
field, switching involving absorbed photons is impractical at a nanoscale). 
Indeed, the applied electric field, on the order of a typical breakdown field 
Eb < lO^V/cm, is much smaller than a typical atomic field ~ lO^V/cm, char- 
acteristic of the energy barriers. Small barrier would be a subject for sporadic 
thermal switching, whereas a larger barrier ~ 1 — 2eV would be impossible to 
overcome with the applied field. One may only change the relative energy of 
the minima by external field and, therefore, redistribute the molecules statisti- 
cally slightly inequivalently between the two states. An intrinsic disadvantage 
of the conformational mechanism, involving motion of ionic group, exceeding 
the electron mass by many orders of magnitude, is a slow switching speed 
(~kHz). In case of supramolecular complexes like rotaxanes and catenanes 
[55] there are two entangled parts which can change mutual positions as a 
result of redox reactions (in solution). Thus, for the rotaxane-based memory 
devices a slow switching speed of ~ 10~^ seconds was reported. 

As a possible conformational E-tie\d addressable molecular switch, we have 
considered a bistable molecule with -CONH2 dipole group [73]. The barrier 
height is Ei, = O.lSeV. Interaction with an external electric field changes the 
energy of the minima, but estimated switching field is huge, ~ 0.5V/A. At 
non-zero temperatures, temperature fiuctuations might result in statistical 
dipole fiipping at lower fields. The I-V curve shows hysteresis in the 3 to 4 
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(b) Current hysteresis 
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Fig. 10. Schematic ol the molecular quantum dot with central conjugated unit sep- 
arated from the electrodes by wide-band insulating molecular groups. First electron 
tunnels into the dot and occupies an empty (degenerate) state there. If the inter- 
action between the first and second incoming electron is repulsive, (7 > 0, then the 
dot will be in a Coulomb blockade regime (a) . If the electrons on the dot effectively 
attract each other, [7 < 0, the system will show current hysteresis (b). 



Volts window for two possible conformations. One can estimate the thermal 
stability of the state as 58 ps at room temperature, and 33 ms at 77 K. 

We explored a possibility for a fast molecular switching where switching 
is due to strong correlation effects on the molecule itself, so-called molecular 
quantum dot (MQD). The molecular quantum dot consists of a central con- 
jugated unit (containing half-occupied, and, therefore, extended tt— orbitals). 
Fig. [ini Frequently, those are formed from the p-states on carbon atoms, 
which are not saturated (i.e. they do not share electrons with other atoms 
forming strong tr— bonds, with typical bonding-antibonding energy difference 
about IRy). Since the tt— orbitals are half-occupied, they form the HOMO- 
LUMO states. The size of the HOMO-LUMO gap is then directly related 
to the size of the conjugated region d, Fig. [TUl by a standard estimate 
■Ehomo-lumo ~ fi^ /md^ '-^ 2 — b eV. It is worth noting that in the con- 
jugated linear polymers like polyacetylene (— C — C)„ the spread of the 
TT— electron would he d — oo and the expected i^HOMO-LUMO = 0. However, 
such a one-dimensional metal is impossible, Peierls distortion (C=C bond 
length dimerization) sets in and opens up a gap of about ^ 1.5eV at the 
Fermi level|17|. In a molecular quantum dot the central conjugated part is 
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separated from electrodes by insulating groups with saturated cr— bonds, like 
e.g. the alkane chains, Fig. 3. Now, there are two main possibilities for carrier 
transport through the molQD. If the length of at least one of the insulat- 
ing groups ii(2) is not very large (a conductance Gi(2) is not much smaller 
than the conductance quantum Go = 2e^//i), then the transport through the 
MQD will proceed by resonant tunneling processes. If, on the other hand, 
both groups are such that the tunnel conductance Gi(2) ^ Go, the charge 
on the dot will be quantized. Then we will have another two possibilities: (i) 
the interaction of the extra carriers on the dot is repulsive U > 0, and we 
have a Coulomb blockade |75] . or (ii) the effective interaction is attractive, 
U < 0, then we would obtain the current hysteresis and switching (see be- 
low) . Coulomb blockade in molecular quantum dots has been demonstrated in 
Refs. [76]. In these works, and in Ref. [77], the three-terminal active molecular 
devices have been fabricated and successfully tested. 

Much faster switching compared to the conformational one, described in 
the previous section, may be caused by coupling to the vibrational degrees 
of freedom, if the vibron-mediated attraction between two carriers on the 
molecule is stronger than their direct Coulomb repulsion [52], Fig. [TOb. The 
attractive energy is the difference of two large interactions, the Coulomb re- 
pulsion and the phonon mediated attraction, on the order of leV each, hence 
\U\ - O.leV. 

5.2 Polaron mechanism of carrier attraction and switching 

Although the correlated electron transport through mesoscopic systems with 
repulsive electron-electron interactions received considerable attention in the 
past, and continues to be the focus of current studies, much less has been 
known about a role of electron-phonon correlations in "molecular quantum 
dots" (MQD). Some while ago we have proposed a negative— L/ Hubbard 
model of a d-fold degenerate quantum dot [78] and a polaron model of resonant 
tunneling through a molecule with degenerate level [52:. We found that the 
attractive electron correlations caused by any interaction within the molecule 
could lead to a molecular switching effect where I-V characteristics have two 
branches with high and low current at the same bias voltage. This prediction 
has been confirmed and extended further in our theory of correlated transport 
through degenerate MQDs with a full account of both the Coulomb repulsion 
and realistic electron-phonon (e-ph) interactions [5^. We have shown that 
while the phonon side-bands significantly modify the shape of hysteretic I-V 
curves in comparison with the negative-f/ Hubbard model, switching remains 
robust. It shows up when the effective interaction of polarons is attractive and 
the state of the dot is multiply degenerate, d> 2. 

First, we shall describe a simplest model of a single atomic level coupled 
with a single one-dimensional oscillator using the first quantization represen- 
tation for its displacement x [79] . 
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Fig. 11. Two localized electrons at sites 1 and 2 with opposite spins shift the 
equilibrium position of the ion at site 3. As a result, the two electrons attract each 
other. 



H = eon + f.n-—— + —. (9) 

Here M and k are the oscillator mass and the spring constant, / is the inter- 
action force, and Ti — c = ks = ^- This Hamiltonian is readily diagonahzed 
with the exact displacement transformation of the vibration coordinate x, 

x^y-nf/k, (10) 

to the transformed Hamiltonian without electron-phonon coupling, 

u - I ky-" 

H — en H , (H) 

2M dy^ 2 ' ^ ' 

e^eo- Ep, (12) 

where we used — fi because of the Fermi-Dirac statistics. It describes a 
small polaron at the atomic level Eq shifted down by the polaron level shift 
Ep — p /2k^ and entirely decoupled from ion vibrations. The ion vibrates near 
a new equilibrium position, shifted by //fc, with the "old" frequency (A:/M)^/^. 
As a result of the local ion deformation, the total energy of the whole system 
decreases by Ep since a decrease of the electron energy by —2Ep overruns an 
increase of the deformation energy Ep. Note that the authors of Ref. [80] made 
an illegitimate replacement of the square of the occupation number operator 
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n = CqCq in Eq. (jlip by its "mean-field" expression — uqu which contains 
the average population of a single molecular level, uq, in disagreement with 
the exact identity, = fi. This leads to a spurious self-interaction of a single 
polaron with itself [i.e. the term e = eo ~ n^Ep instead of Eq. P^ ]. and a 
resulting non-existent nonlinearity in the rate equation. 

Lattice deformation also strongly affects the interaction between electrons. 
When a short-range deformation potential and molecular e-ph interactions are 
taken into account together with the long-range Frohlich interaction, they can 
overcome the Coulomb repulsion. The resulting interaction becomes attractive 
at a short distance comparable to a lattice constant. The origin of the attrac- 
tive force between two small polarons can be readily understood from a similar 
Holstein-like toy model as above [81], but with two electrons on neighboring 
sites 1,2 interacting with an ion 3 between them. Fig. [TT] For generality, we 
now assume that the ion is a three-dimensional oscillator described by a dis- 
placement vector u, rather than by a single-component displacement x as in 
Eq.(l). 

The vibration part of the Hamiltonian in the model is 

Electron potential energies due to the Coulomb interaction with the ion are 
about 

^1,2 =14)-uVfl,,,,yo(i?l,2), (14) 

where R.i(2) is the vector connecting ion site 3 with electron 1(2). Hence, the 
Hamiltonian of the model is given by 



H = Ea{ni + ns) + u • (fifii + hn2) - T^Iti^ + (15) 



where fi_2 — -^6^61.2/0^ is the Coulomb force, and n\^2 are occupation number 
operators at every site. This Hamiltonian is also readily diagonalizcd by the 
same displacement transformation of the vibronic coordinate u as above, 

u = v-(firii+f2n2)/fc. (16) 

The transformed Hamiltonian has no electron-phonon coupling, 

19^ fcv^ 

H = {eo- Ep){ni + fia) + Vphnin2 ~ + — ' (I'') 

and it describes two small polarons at their atomic levels shifted by the polaron 
level shift Ep = /^2/2^, which are entirely decoupled from ion vibrations. As 
a result, the lattice deformation caused by two electrons leads to an effec- 
tive interaction between them, Vph, which should be added to their Coulomb 
repulsion, Vc, 

Vph = -fi • i2/k. (18) 
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When Vph is negative and larger by magnitude than the positive Vc, the re- 
sulting interaction becomes attractive. That is Vph rather than Ep, which is 
responsible for the hysteretic behavior of MQDs, as discussed below. 



5.3 Exact solution 



The procedure, which fully accounts for all correlations in MQD is as follows, 
see Ref. [52]. The molecular Hamiltonian includes the Coulomb repulsion, U'^ , 
and the electron-vibron interaction as 



H = e n +- U^n'fi flu' 

+ II «.^9(7m9C«9 + ^-c-) + '^M^i + 1/2)- (19) 



Here dq annihilates phonons, ujq is the phonon (vibron) frequency, and 
are the e-ph coupling constant (g enumerates the vibron modes). This Hamil- 
tonian conserves the occupation numbers of molecular states ri^. 
One can apply the canonical unitary transformation e'^, with 

S ^-^Uf, {jf,qdq - H.C.) 



integrating phonons out. The electron and phonon operators are transformed 
as 



= exp - 
\ 9 



= c^^p, = exp I 2^ iM^q - H.C. ) (20) 

and 

dq ^ dq -^h^,-flq, (21) 
M 

respectively. This Lang-Firsov transformation shifts ions to new equilibrium 
positions with no effect on the phonon frequencies. The diagonalization is 
exact: 

H = I]e>M + II^«Krf« + 1/2) + \Y. '^mm'^mV' (22) 

i q fi^fj.' 

where 

f^MM'^C^M'M'-2ETM97M'9^9, (23) 

q 

is the renormalized interaction of polarons comprising their interaction via 
molecular deformations (vibrons) and the original Coulomb repulsion, J7^/. 
The molecular energy levels are shifted by the polaron level-shift due to the 
deformation created by the polaron. 
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If we assume that the couphng to the leads is weak, so that the level width 
r <^ \U\, we can find the current from 1821 



/oo 
duj[h{oj)- hicu)] p{uj), 
-OO 

pH = -i^ImG», 



(25) 
(26) 



where is a complete set of one-particle molecular states, fi(^2){^) - 
^exp '^+'^^<^^/^ Fermi function. Here Iq ~ qF, p{uj) is the molec- 

ular DOS, Gj}{uj) is the Fourier transform of the Green's function Gjj{t) - 
—iO{t) ^{c^t(i), cji}) ,{•••,•••} is the anticomniutator, Cp(t) = e*^*c^e~*^* 
9{t) = 1 for t > and zero otherwise. We calculate p{uj) exactly for the Hamil- 
tonian (|22p. which includes both the Coulomb U'~^ and e-ph interactions. 
The retarded GF becomes 



G^^^t)^~^e{t)[{c,{t)cl){X,{t)Xl) 
+ {elegit)) {XlX^it))]. 

The phonon correlator is simply 

{X,{t)Xl)^e^^Y.- 



(27) 



q sinh^ 



cos I Lut + i 



cosh ^— 2- 
2 



(28) 



where the inverse temperature (3 — 1/T, and (^Xj^X^{t)^ — (X^(t)X^j)* . 
The remaining GFs (c^{t)c^^'^, are found from the equations of motion ex- 
actly. Finally, for the simplest case of a coupling to a single mode with the 
characteristic frequency wq and 7^ = 7 we find [S2] 



G«(c.)=Z^a(n)^/aC) 



r=0 



e 2 



1 



uj — rU — IljJq + iS uj — rU + Iloq + id 



(l-<5io)e 

1 — n 



rU + lujQ + i8 Lo — rU — Iujq + iS 



(29) 



where 
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(30) 



q 

is the familiar polaron narrowing factor, the degeneracy factor 

Crin)^^j^^^n^l-nr-^-^, (31) 

C = ItP/ sinh^^l^, /; (^) the modified Bessel function, and Sik the Kroneker 
symbol. The important feature of the DOS, Eq. ([^5]) . is its nonlinear depen- 
dence on the occupation number n. It contains full information about all 
possible correlation and inelastic effects in transport, in particular, all the 
phonon sidebands. 

To simplify our discussion, we assume that the Coulomb integrals do not 
depend on the orbital index, i.e. U^^/ = U , and consider a coupling to a single 
vibrational mode, ujq — loq. Applying the same transformation to the retarded 
Green's function, one obtains the exact spectral function [52] for ad— fold 
degenerate MQD (i.e. the density of molecular states, DOS) as 

d—l oo 
r=0 1=0 

„/3^oi/2 _ ^rU- lujo) + nS{u -rU + Icua)] 

+ (1 - 6io)e~^'^°^/^[nS{uj - rU - Iluq) 
+ (1 - n)S{uj -rU + lujo)] 



(32) 



where Z — exp — |7pcoth^^|^ , is the above polaron narrowing factor, ^ — 

|7p/ sinh(/3Li;o/2), /3 = ^/T, and Sik the Kroneker symbol. 

The important feature of DOS, Eq. ([5^ . is its nonlinear dependence on 
the average electronic population n — (c)^c^^ , which leads to the switching, 
hysteresis, and other nonlinear effects in I-V characteristics for d > 2 [52 . It 
appears due to correlations between different electronic states via the corre- 
lation coefficients Cr{n). There is no nonlinearity if the dot is nondegenerate, 
d = 1, since Co(n) = 1. In this simple case the DOS, Eq. ([5^ . is a linear 
function of the average population that can be found as a textbook example 
of an exactly solvable problems [83] . 

In the present case of MQD weakly coupled with leads, one can apply 
the Fermi-Dirac golden rule to obtain an equation for n. Equating incoming 
and outgoing numbers of electrons in MQD per unit time we obtain the self- 
consistent equation for the level occupation n as 
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/CO 
du;{rih{iu) + r2f2{uj)}p{iu) 
-oo 

du;{r,[l - /iH] +r2[i - /2H]}pM, (33) 



where A(2) are the transition rates from left (right) leads to MQD, and p{llj) 
is found from Eqs. and (|25p . For d = 1,2 the kinetic equation for n is 
linear, and the switching is absent. Switching appears for d > 3, when the 
kinetic equation becomes non-linear. Taking into account that p{^) = 
Eq. ([55)1 for the symmetric leads, Fi = I2, reduces to 

2nd^ J dLop{u){h + h), (34) 

which automatically satisfies < n < 1. Explicitly, the self-consistent equation 
for the occupation number is 

1 "^"^ 

^ -Y^Cr{n)[na+ + [l - n)b+], (35) 



n 

2 

r=0 



a+ = Z ^ (0 ( [h{rU ~ Icoo) + h{rU - lu:^)] 



where 



+ (1 - 6iQ)e-— [h{rU + luj^) + ^(rf/ + lujo)]J , (36) 

ht = zY^hiOie^ [/i {rU + lu;o) + /a (rf/ + ^c^o)] 
1=0 ^ 

+ (1 - (5jo)e-^ [/i(r[/ - Zcjo) + /2(r[/ - Iluq)]^ . (37) 
The current is expressed as 



where 



a- =ZJ2 Ii (0 ( e ^ [/i (^C^ - ^^o) - /2 {rU - IcJo)] 



+ (1 - <5jo)e-^ [fi{rU + lujo) ~ f2{rU + lujo)]^ , (39) 

OC / 

6- = Z ^ (0 ( [/i (rC/ -f ^t^o) - /2 (rC/ + Zc^o)] 

(1 - Sio)e-^ [fi{rU ~ Icuo) - f2{rU - Iljo)]) . (40) 
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Fig. 12. Current-voltage characteristic of the nondegenerate (d — 1) MQD at T = 0, 
oJo/A = 0.2, and j'^ — 11/13. There is the phonon ladder in I-V , but no hysteresis. 



5.4 Absence of switching of single- or double-degenerate MQD 

If the transition rates from electrodes to MQD are small, <C wq, the rate 
equation for n and the current, I{V) are readily obtained by using the exact 
molecular DOS, Eq. ([5^ and the Fermi-Dirac Golden rule. In particular, for 
the nondegenerate MQD and T = OK the result is 
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The general expressions for the coefhcients are given at arbitrary temperatures 
in Ref. [52j . and they are especially simple in low-temperature limit: 
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Fig. 13. Current-voltage characteristic of two-fold degenerate MQDs (d — 2) does 
not show hysteretic behavior. The parameters are the same as in Fig. 1121 Larger 
number of elementary processes for conductance compared to the nondegenerate 
case of d = 1 generates more steps in the phonon ladder in comparison with Fig. 1121 
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j — I/Iq, Iq — edr, A is the position of the MQD level with respect to the 
Fermi level at V = 0, and 6{x) = 1 if x > and zero otherwise. The current is 
single valued, Fig. [121 with the familiar steps due to phonon-side bands. The 
double-degenerate level provides more elementary processes for conductance 
reflected in larger number of steps on phonon ladder compared to rf = 2 case, 
Fig.IIl 
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On the contrary, the mean-field approximation (MFA) leads to the op- 
posite conclusion. Galperin et al. [5D] have replaced the occupation number 
operator fi in the e-ph interaction by the average population ng [Eq. (2) of 
Ref. [5D]] and found the average steady-state vibronic displacement {d + d)) 
proportional to no (this is an explicit neglect of all quantum fluctuations on 
the dot accounted for in the exact solution). Then, replacing the displacement 
operator d -f in the bare Hamiltonian, Eq. (11), by its average, Ref. [80] . 
they obtained a new molecular level, eo = £o — 2£reorg?T-o shifted linearly with 
the average population of the level. This is in stark disagreement with the 
conventional constant polaronic level shift, Eq. (|12l24p {ereorg is |7pWo in 
our notations). The MFA spectral function turned out to be highly nonlin- 
ear as a function of the population, e.g. for the weak-coupling with the leads 
p{ijj) — 5{lo — £o — '^£reorgno), scc Eq. (17) in Ref. [80j. As a result, the au- 
thors of Ref. [50] have found multiple solutions for the steady-state population, 
Eq. (15) and Fig. 1, and switching. Fig. 4 of Ref. [5D], which actually do not 
exist being an artefact of the model. 

In the case of a double-degenerate MQD, d — 2, there are two terms, which 
contribute to the sum over r, with Co(n) = 1 — n and Ci(n) = n. The rate 
equation becomes a quadratic one [52] . Nevertheless there is only one physical 
root for any temperature and voltage, and the current is also single- valued. 
Fig. 3. 

Note that the mean-field solution by Galperin et al. [80] applies at any ra- 
tio r /u)Q, including the limit of interest to us, F <C wo- where their transition 
between the states with tiq = and 1 only sharpens, but none of the results 
change. Therefore, MFA predicts a current bistability in the system where it 
does not exist at d = 1. Ref. [80] plots the results for F > wq, ~ 0.1 — 0.3 
eV, which corresponds to molecular bridges with a resistance of about a few 
100 kl2. Such model "molecules" are rather "metallic" in their conductance 
and could hardly show any bistability at all because carriers do not have time 
to interact with vibrons on the molecule. Indeed, taking into account the cou- 
pling with the leads beyond the second order and the coupling between the 
molecular and bath phonons could hardly provide any non-linearity because 
these couplings do not depend on the electron population. This rather obvi- 
ous conclusion for molecules strongly coupled to the electrodes can be reached 
in many ways, see e.g. a derivation in Refs. [HU |85]. While Refs. [HU [85] do 
talk about telegraph current noise in the model, there is no hysteresis in the 
adiabatic regime, F ^ wq either. This result certainly has nothing to do 
with our mechanism of switching [52j that applies to molecular quantum dots 
{F <C Wo) with d> 2. Such regime has not been studied in Refs. [84 ] [85 ] [86] . 
which have applied the adiabatic approximation, as being "too challenging 
problem". Nevertheless, Mitra et al. [5S] have misrepresented our formalism 
[52] claiming that it "lacks of renormanlization of the dot-lead coupling" (due 
to electron- vibron interaction), or "treats it in an average manner". In fact, 
the formalism [52j is exact, fully taking into account the polaronic renor- 
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malization, phonon-side bands and polaron-polaron correlations in the exact 
molecular DOS, Eq. ([5^ . 

As a matter of fact, most of the molecules are very resistive, so the actual 
molecular quantum dots are in the regime we study, see Ref . [87| . For exam- 
ple, the resistance of fully conjugated three-phenyl ring Tour-Reed molecules 
chemically bonded to metallic Au electrodes [33] exceeds lGf2. Therefore, 
most of the molecules of interest to us are in the regime that we discussed, 
not that of Refs.[84l[85]. 

5.5 Nonlinear rate equation and switching 

Now, consider the case d = 4, where the rate equation is non-linear, to see if it 
produces multiple physical solutions. For instance, in the limit I7I <C 1, T = 0, 
where we have br ^ Ur, Z = 1, the remaining interaction is U = < 0, we 
recover the negative-?/ model [75|, and the kinetic equation for d = 4 is 

2n = 1 - (1 - nf (45) 

in the voltage range Z\ — |[/| < eV/2 < A. This equation has an additional 
nontrivial physical solution n = {3—^/5)/2 = 0.38. The current is simplified as 
I/Io = 2n. The current- voltage characteristics will show a hysteretic behavior 
in this case for d — 4. When the voltage increases from zero, 4-fold degenerate 
MQD remains in a low-current state until the threshold eV2/2 = Z\ is reached. 
Remarkably, when the voltage decreases from the value above the threshold 
V2, the molecule remains in the high-current state down to the voltage eVi/2 — 
A-\U\ well below the threshold ¥2- 

In fact, there is a hysteresis of current at all values of the electron-phonon 
constant 7, e.g. 7^ = 11/13 (selected in order to avoid an accidental com- 
mensurability of ujq and U), Fig. [141 Indeed, the exact equation for average 
occupation of the dot reads 

2n = (1 - n f[na+ + (1 - n)b+] 

+3n{l - nf[na+ + (1 - n)b+] 
+3n^{l - n)[na+ + (1 - n)6^] 

W[na^ + il-n)bj]. (46) 

We solved this nonlinear equation for the case llIq/A — 0.2, U'~'' = 0, so 
that the attraction between electrons is U = —2j'^u!o — —0.4 (all energies 
are in units of A), and found an additional stable solution for an average 
occupation number (and current) that results in a hysteresis curve. Fig. 1141 
The bistability region shrinks down with temperature, and the hysteresis loop 
practically closes at T/A = 0.01. As we see from Eqs. (|36| . (jST j) . the electron 
levels with phonon sidebands Z\±/a;o, A + Uztlu!o^ A + 2U±Iluo, A + 3U±Iujo 
with / = 0, 1, ... contribute to electron transport with different weights, and 
this creates a complex picture of steps on the I-V curve. Fig. [TH 
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Fig. 14. The I-V curves for tunneling through the molecular quantum dot, Fig. llOb 
with the electron-vibron coupling constant 7^ — 11/13. 
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Note that switching required a degenerate MQD (c? > 2) and the weak 
couphng to the electrodes, F <^ ljq. Different from the non-degenerate dot, 
the rate equation for a multi-degenerate dot, d > 2, weakly coupled to the 
leads has multiple physical roots in a certain voltage range and a hysteretic 
behavior due to correlations between different electronic states of MQD [52] . 

Summarizing this Section, we have calculated the I-V characteristics of the 
nondegenerate and two-fold degenerate MQDs showing no hysteretic behav- 
ior, and conclude that mean field approximation [80] leads to a non-existent 
hysteresis in a model that was solved exactly in Ref. [S2]. Different from the 
non-degenerate and two-fold degenerate dots, the rate equation for a multi- 
degenerate dot, d> 2, weakly coupled to the leads, has multiple physical roots 
in a certain voltage range showing hysteretic behavior due to correlations be- 
tween different electronic states of MQD [52] . Our conclusions are important 
for searching of the current-controlled polaronic molecular switches. Inciden- 
tally, Ceo molecules have the degeneracy d = 6 of the lowest unoccupied 
level, which makes them one of the most promising candidate systems, if the 
weak-coupling with leads is secured. 

6 Role of defects in molecular transport 

Interesting behavior of electron transport in molecular systems, as described 
above, refers to ideal systems without imperfection in ordering and com- 
position. In reality, one expects that there will be a considerable disorder 
and defects in organic molecular films. As mentioned above, the conduction 
through absorbed [54 and Langmuir-Blodgett [55] monolayers of fatty acids 
(CIl2)„ was associated with defects. An absence of tunneling through self- 
assembled monolayers of C12-C18 (inferred from an absence of thickness de- 
pendence at room temperature) has been reported by Boulas et al. [53]. On 
the other hand, the tunneling in alkanethiol SAMs was reported in [88] I89j . 
with an exponential dependence of monolayer resistance on the chain length 
L, Ra cc ex.p(f3crL), and no temperature dependence of the conductance in 
C8-C16 molecules was observed over the temperatures T = 80 — 300K [89] . 

The electrons in alkane molecules are tightly bound to the C atoms by 
£7— bonds, and the band gap (between the highest occupied molecular orbital, 
HOMO, and lowest unoccupied molecular orbital, LUMO) is large, ~ 9 — lOeV 
[53] . In conjugated systems with tt— electrons the molecular orbitals are ex- 
tended, and the HOMO-LUMO gap is correspondingly smaller, as in e.g. 
polythiophenes, where the resistance was also found to scale exponentially 
with the length of the chain, Rt^ cx exp(/3^L), with /J^r — 0.35A^^ instead of 
Per = 1.08A~^[88^. In stark contrast with the temperature-independent tun- 
neling results for SAMs [89], recent extensive studies of electron transport 
through 2.8 nm thick eicosanoic acid (C20) LB monolayers at temperatures 
2K-300K have established that the current is practically temperature inde- 
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Fig. 15. Local density of states and transmission as a function of energy for (a) C13 
with Au impurity and (b) C13 with Au impurity and H vacancy (dangling bond). 
Middle sections show closeups of the resonant peaks due to deep defect levels with 
respect to the HOMO and LUMO molecular states. The HOMO-LUMO gap is about 
10 eV. 



pendent below T < 60K, but very strongly temperature dependent at higher 
temperatures T = 60 - 300K[90]. 

A large amount of effort went into characterizing the organic thin films 
and possible defects there [TSl [HIl ■ It has been found that the electrode 
material, like gold, gets into the body of the film, leading to the possibility 
of metal ions existing in the film as single impurities and clusters. Electronic 
states on these impurity ions are available for the resonant tunneling of car- 
riers in very thin films (or hopping in thicker films, a crossover between the 
regimes depending on the thickness). Depending on the density of the impu- 
rity states, with increasing film thickness the tunneling will be assisted by 
impurity "chains" , with an increasing number of equidistant impurities [93j . 
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One-impurity channels produce steps on the I-V curve but no temperature de- 
pendence, whereas the inelastic tunneling through pairs of impurities at low 
temperatures defines the temperature dependence of the film conductance, 
G(T) oc T''/^, and the voltage dependence of current I{V) cx This 
behavior has been predicted theoretically and observed experimentally for 
tunneling through amorphous Si [95] and Al2 03[96l. Due to the inevitable 
disorder in a "soft" matrix, the resonant states on different impurities within 
a "channel" will be randomly moving in and out of resonance, creating meso- 
scopic fluctuations of the I-V curve. The tunneling may be accompanied by 
interaction with vibrons on the molecule, causing step-like features on the I-V 
curve [ZZlEl]. 

During processing, especially top electrode deposition, small clusters of the 
electrode material may form in the organic film, causing Coulomb blockade, 
which also can show up as steps on the TV curve. It has long been known that 
a strong applied field can cause localized damage to thin films, presumably 
due to electromigration and the formation of conducting filaments [56] . The 
damaged area was about 30nm in diameter in 40-160 monolayer thick LB 
films [56] (a) and 5-10/im in diameter in films 500-5000A thick, and showed 
switching behavior under external bias voltage cycling [56] (b). As discussed 
above, recent spatial mapping of a conductance in LB monolayers of fatic 
acids with the use of conducting AFM has revealed damage areas 30-lOOnm in 
diameter, frequently appearing in samples after a "soft" electrical breakdown, 
which is sometimes accompanied by a strong temperature dependence of the 
conductance through the film[58]. 

A crossover from tunneling at low temperatures to an activation-like de- 
pendence at higher temperatures is expected for electron transport through 
organic molecular films. There are recent reports about such a crossover in in- 
dividual molecules like the 2nm long Tour wire with a small activation energy 
Ea ~ ISOmeV [97]. Very small activation energies on the order of 10-100 meV 
have been observed in polythiophene monolayers [HH] • Our present results sug- 
gest that this may be a result of interplay between the drastic renormalization 
of the electronic structure of the molecule in contact with electrodes, and dis- 
order in the film (Fig. [HI right inset). We report the ab-initio calculations 
of point-defect assisted tunneling through alkanedithiols S(CH2)„S and thio- 
phene T3 (three rings SC4) self-assembled on gold electrodes. The length of 
the alkane chain was in the range n = 9 — 15. 

We have studied single and double defects in the film: (i) single Au im- 
purity. Figs. 10a, 11a, (ii) Au impurity and H vacancy (dangling bond) on 
the chain. Figs. fTBb . fTSh . (iii) a pair of Au impurities. Fig. [T5b . (iv) Au 
and a "kink" on the chain (one C=C bond instead of a C-C bond). Single 
defect states result in steps on the associated I-V curve, whereas molecules 
in the presence of two defects generally exhibit a negative differential re- 
sistance (NDR). Both types of behavior are generic and may be relevant 
to some observed unusual transport characteristics of SAMs and LB films 
|54l [55l [SHI [90l [97l [98] . We have used an ab-initio approach that combines 
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the Keldysh non-equilibrium Green's function (NEGF) method with self- 
consistent pseudopotential-based real space density functional theory (DFT) 
for systems with open boundary conditions provided by semi-infinite elec- 
trodes under external bias voltage [51] I^Hj ■ AH present structures have been 
relaxed with the Gaussian98 code prior to transport calculations [99]. The con- 
ductance of the system at a given energy is found from Eq. ([8]) and the current 
from Eq. (j4|). 

The equilibrium position of an Au impurity is about SA away from the 
alkane chain, which is a typical Van-der-Waals distance. As the density maps 
show (Fig. I16[) . there is an appreciable hybridization between the s- and d- 
states of Au and the sp- states of the carbohydrate chain. Furthermore, the 
Au"*" ion produces a Coulomb center trapping a 6s electron state at an energy 
Ci — — 0.35eV with respect to the Fermi level, almost in the middle of the 
HOMO-LUMO ~ lOeV gap in Cn. The tunneling evanescent resonant state 
is a superposition of the HOMO and LUMO molecular orbitals. Those or- 
bitals have a very complex spatial structure, reflected in an asymmetric line 
shape for the transmission. Since the impurity levels are very deep, they may 
be understood within the model of "short-range impurity potential" jlOOj . 
Indeed, the impurity wave function outside of the narrow well can be fairly 
approximated as 

= (47) 

where k is the inverse radius of the state, h^K,^/2m* = Ei, where Ei — A — ti is 
the depth of the impurity level with respect to the LUMO, and Z\=LUMO— F 
is the distance between the LUMO and the Fermi level F of gold and, con- 
sequently, the radius of the impurity state 1/k is small. The energy distance 
A K, 4.8eV in alkane chains (CH2)„[53] (« 5eV from DFT calculations), and 
m* ^ 0.4 the effective tunneling mass in alkanes[89]. For one impurity in a 
rectangular tunnel barrier [100] we obtain the Breit-Wigner form of transmis- 
sion T{E, V), as before, Eq. ([5]). Using the model with the impurity state wave 
function (|Tf)) . we may estimate for an Au impurity in C13 [L = 10. 9A) the 
width Fl = Ffi = 1.2 x 10~^, which is within an order of magnitude compared 
with the calculated value 1.85 x 10~^eV. The transmission is maximal and 
equals unity when E — ei and F^ — F^i, which corresponds to a symmetrical 
position for the impurity with respect to the electrodes. 

The electronic structure of the alkane backbone, through which the elec- 
tron tunnels to an electrode, shows up in the asymmetric lineshape, which is 
substantially non-Lorentzian, Fig. 1151 The current remains small until the bias 
has aligned the impurity level with the Fermi level of the electrodes, resulting 
in a step in the current, h ~ ^FQe~'^^ fFig.[T5b). This step can be observed 
only when the impurity level is not very far from the Fermi level F, such that 
biasing the contact can produce alignment before a breakdown of the device 
may occur. The most interesting situations that we have found relate to the 
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Fig. 16. Current-voltage characteristics of an alkane chain C13 with (a) single 
An impurity (6s-state), (b) two Au impurities (5d and 6s-states on left and right 
ions, respectively), and (c) Au impurity and H vacancy (dangling bond). Double 
defects produce the negative differential resistance peaks (b) and (c). Inset shows 
the density of states, transmission, and stick model for polythiophene T3. There is 
significant transmission at the Fermi level, suggesting an ohmic I-V characteristic 
for T3 connected to gold electrodes. Disorder in the film may localize states close to 
the Fermi level (schematically marked by arrow), which may assist in hole hopping 
transport with an apparently very low activation energy (0.01-O.leV), as is observed. 
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pairs of point defects in the film. If the concentration of defects is c ^ 1, the 
relative number of configurations with pairs of impurities will be very small, 
(X c^. However, they give an exponentially larger contribution to the current. 
Indeed, the optimal position of two impurities is symmetrical, a distance L/2 
apart, with current I2 oc e^"^/^. The conductance of a two-impurity chain is 

m 

9i2(^) - —r-— -—-TZ, — ::r72- (-48) 



\iE - ei + irL)iE - 62 + iFii) - f- 



121 



For a pair of impurities with slightly differing energies ti2 = 2(£'i+i?2)e^'"'^^/Kri2, 
where ri2 is the distance between them. The interpretation of the two-impurity 
channel conductance (|48p is fairly straightforward: if there were no coupling 
to the electrodes, i.e. Fl = Fjf = 0, the poles of gi2 would coincide with the 
bonding and antibonding levels of the two-impurity "molecule" . The coupling 
to the electrodes gives them a finite width and produces, generally, two peaks 
in conductance, whose relative positions in energy change with the bias. The 
same consideration is valid for longer chains too, and gives an intuitive picture 
of the formation of the impurity "band" of states. The maximal conductivity 
gi2 — lirfi occurs when ei = £2, r^ri^. = ^12 = ^21 where /2 is the width of 
the two-impurity resonance, and it corresponds to the symmetrical position of 
the impurities along the normal to the contacts separated by a distance equal 
to half of the molecule length, ri2 = L/2. The important property of the 
two- impurity case is that it produces negative differential resistance (NDR). 
Indeed, under external bias voltage the impurity levels shift as 

Ci = do + qVzi/L, (49) 

where Zi are the positions of the impurity atoms counted from the center of 
the molecule. Due to disorder in the film, under bias voltage the levels will 
be moving in and out of resonance, thus producing NDR peaks on the I-V 
curve. The most pronounced negative differential resistance is presented by a 
gold impurity next to a Cn chain with an H- vacancy on one site. Fig. 115b (the 
defect corresponds to a dangling bond). The defects result in two resonant 
peaks in transmission. Surprisingly, the H vacancy (dangling bond) has an 
energy very close to the electrode Fermi level F, with Ci = — O.leV (Fig. [TSb. 
right peak). The relative positions of the resonant peaks move with an external 
bias and cross at 1.2V, producing a pronounced NDR peak in the I-V curve, 
Fig. \Wb . No NDR peak is seen in the case of an Au impurity and a kink 
C=C on the chain because the energy of the kink level is far from that of 
the Au 6s impurity level. The calculated values of the peak current through 
the molecules were large: Ip « 90 nA/molecule for an Au impurity with H 
vacancy, and ~ 5 nA/molecule for double Au impurities. 

We have observed a new mechanism for the NDR peak in a situation with 
two Au impurities in the film. Namely, Au ions produce two sets of deep 
impurity levels in Cn films, one stemming from the 6s orbital, another from 
the 5d shell, as clearly seen in Fig. [T6b (inset). The 5d-states are separated 
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in energy from 6s, so that now the tunnehng through s-d pairs of states is 
allowed in addition to s-s tunneling. Since the 5d-states are at a lower energy 
than the s-state, the d- and s-states on different Au ions will be aligned at 
a certain bias. Due to the different angular character of those orbitals, the 
tunneling between the s-state on the first impurity and a d-state on another 
impurity will be described by the hopping integral analogous to the Slater- 
Koster sda integral. The peak current in that case is smaller than for the pair 
Au-H vacancy, where the overlap is of ssa type (cf. Figs. fTBb.c). 

Thiophene molecules behave very differently since the tt— states there are 
conjugated and, consequently, the HOMO-LUMO gap is much narrower, just 
below 2 eV. The tail of the HOMO state in the T3 molecule (with three 
rings) has a significant presence near the electrode Fermi level, resulting in 
a practically "metallic" density of states and hence ohmic I-V characteristic. 
This behavior is quite robust and is in apparent disagreement with experiment, 
where tunneling has been observed[55j. However, in actual thiophene devices 
the contact between the molecule and electrodes is obviously very poor, and 
it may lead to unusual current paths and temperature dependence [98]. 

We have presented the first parameter-free DFT calculations of a class 
of organic molecular chains incorporating single or double point defects. The 
results suggest that the present generic defects produce deep impurity levels in 
the film and cause a resonant tunneling of electrons through the film, strongly 
dependent on the type of defects. Thus, a missing hydrogen produces a level 
(dangling bond) with an energy very close to the Fermi level of the gold 
electrodes F. In the case of a single impurity, it produces steps on the I-V 
curve when one electrode's Fermi level aligns with the impurity level under 
a certain bias voltage. The two-defect case is much richer, since in this case 
we generally see a formation of the negative-differential resistance peaks. We 
found that the Au atom together with the hydrogen vacancy (dangling bond) 
produces the most pronounced NDR peak at a bias of 1.2V in C13. Other 
pairs of defects do not produce such spectacular NDR peaks. A short range 
impurity potential model reproduces the data very well, although the actual 
lineshape is different. 

There is a remaining question of what may cause the strong tempera- 
ture dependence of conductance in "simple" organic films like [CH2]„. The 
activation-like conductance cx exp{—Ea/T) has been reported with a small 
activation energy Ea ~ 100 — 200meV in alkanes [221 HZ] and even smaller, 
10-100 meV, in polythiophenes [98]. This is much smaller than the value calcu- 
lated here for alkanes and expected from electrical and optical measurements 
on Cn molecules, Ea ^ A ^ 4eV [53], which correspond nicely to the present 
results. In conjugated systems, however, there may be rather natural expla- 
nation of small activation energies. Indeed, the HOMO in T3 polythiophene 
on gold is dramatically broadened, shifted to higher energies and has a con- 
siderable weight at the Fermi level. The upward shift of the HOMO is just a 
consequence of the work function difference between gold and the molecule. In 
the presence of (inevitable) disorder in the film some of the electronic states 
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on the molecules will be localized in the vicinity of Ep. Those states will 
assist the thermally activated hopping of holes within a range of small activa- 
tion energies < 0.1 eV. Similar behavior is expected for Tour wires where 
Ep — HOMO ^ 1 eV[5T](c), if the electrode-molecule contact is poor, as is 
usually the case. 

With regards to carrier hopping in monolayers of saturated molecules, one 
may reasonably expect that in many studied cases the organic films are riddled 
with metallic protrusions (filaments), emerging due to electromigration in a 
very strong electric field, and/or metallic, hydroxyl, etc. inclusions [551 [55] . It 
may result in a much smaller tunneling distance d for the carriers and the 
image charge lowering of the barrier. The image charge lowering of the barrier 
in a gap of width d is AU = g^ln4/(ed), meaning that a decrease of about 
3.5eV may only happen in an unrealistically narrow gap d = 2 — SA in a film 
with dielectric constant e = 2.5, but it will add to the barrier lowering. More 
detailed characterization and theoretical studies along these lines may help to 
resolve this very unusual behavior. We note that such a mechanism cannot 
explain the crossover with temperature from tunneling to hopping reported 
for single molecular measurements, which has to be a property of the device, 
but not a single molecule [97] . 

7 Conclusions 

Studying molecules as possible building blocks for ultradense electronic cir- 
cuits is a fascinating quest that spans of over 30 years. It was inspired decades 
ago by the notion that silicon technology is approaching its limiting feature 
size, estimated at around 1985 to be about 1/im llOl] . More than thirty years 
later and with FET gate lengths getting below lOnm [102j . the same notion 
that silicon needs to be replaced at some point by other technologies floats 
again. We do not know whether alternatives will continue to be steamrollered 
by silicon technology, which is a leading nanotechnology at the moment, but 
the mounting resistance to the famed Moore's law requires to look hard at 
other solutions for power dissipation, leakage current, crosstalk, speed, and 
other very serious problems. There are very interesting developments in study- 
ing electronic transport through molecular films but the mechanisms of some 
observed conductance "switching" and/or nonlinear electric behavior remain 
elusive, and this interesting behavior remains intermittent and not very re- 
producible. Most of the currently observed switching is extrinsic in nature. 
For instance, we have discussed the effect of molecule-electrode contact: the 
tilting of the angle at which the conjugated molecule attaches to the elec- 
trode may dramatically change its conductance, and that probably explains 
extrinsic "telegraph" switching observed in Tour wires [23l [34] and molecule 
reconfigurations may lead to similar phenomena in other systems [ST]. Defects 
in molecular films have also been discussed and may result in spurious peaks in 
I-V curves. We have outlined some designs of the molecules that may demon- 
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strate rectifying behavior, which we call molecular "quantum dots". We have 
shown that at least in some special cases molecular quantum dots may exhibit 
fast (~THz) intrinsic switching. 

The author is grateful to Jeanie Lau, Jason Fitters and Robert Wolkow 
for kind permission to use their data and figures. The work has been partly 
supported by DARPA. 
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